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Abstract 

In this paper we study the effect of normal and shear strains and oxygen vacancies on the structure of 
180° ferroelectric domain walls in PbTiOs. It is known that oxygen vacancies move to the domain walls 
and pin them. Hence, we assume a periodic arrangement of oxygen vacancies on both Pb-centered and Ti- 
centered domain walls in PbTiOs. We use a semi-analytic anharmonic lattice statics method for obtaining 
the relaxed configurations using a shell potential. In agreement with recent ab initio calculations, we observe 
that a Pb-centered domain wall with oxygen vacancies is not stable even under strain. Our semi-analytic 
calculations for PbTiOs show that oxygen vacancies affect the structure of 180° domain walls significantly 
but do not have a considerable effect on the thickness of domain walls; they broaden the domain walls by 
about fifty percent. We also study the effect of normal and shear strains on both perfect and defective 180° 
domain walls. We observe that normal and shear strains affect the structure but do not change the domain 
wall thickness. 

1 Introduction 

Ferroelectric perovskites have been the focus of intense research in recent years because of their potential 
applications in high strain actuators, high density storage devices, etc. [Bhattacharya and Ravichandran, 2003]. 
It is known that macroscopic properties of ferroelcctrics are strongly dependent on domain walls, which are 
extended two-dimensional defects. Any fundamental understanding of ferroelectricity in perovskites requires a 
detailed understanding of domain walls in the nanoscale as these defects are atomically sharp (see Dawber, et al. 
[2005] and references therein). Theoretical studies of domain walls have revealed many of their interesting 
features. From both ab inito calculations [Meyer and Vanderbilt, 2001; Padilla, et al., 1996; Poykko and Chadi, 
1999, 2000] and anaharmonic lattice statics calculations [Yavari, et al., 2007b] it is now known that ferroelectric 
domain walls are atomically sharp. In all these studies, the structure calculations have been done for perfect 
domain walls and free of strain. However, domain walls interact with other types of defects and mainly with 
oxygen vacancies and this may affect the structure and properties of domain walls. Strain may also have a 
significant effect on domain wall structure. 

It is known that the presence of point defects can have important effects on tire properties of perovskites. 
For example, Bujakiewicz-Koronska and Natanzon [2009] show that point defects alter the elastic constants 
of Nai/2Bii/2Ti03 significantly. In this paper, we study the effect of oxygen vacancies and normal and shear 
strains on the structure of 180° ferroelectric domain walls. It is known that oxygen vacancies move to the domain 
walls and pin them and thus we assume a periodic arrangement of oxygen vacancies on both Pb-centered and 
Ti-centered domain walls. We use a semi-analytic anharmonic lattice statics method for obtaining the relaxed 
configurations using a shell potential. In agreement with recent ab initio calculations, we observe that a Pb- 
centered domain wall with oxygen vacancies is not stable even under strain. Our semi-analytic solutions for 
PbTiOa show that oxygen vacancies increase the thickness of the domain wall by about fifty percent. This is 
different from the results of a recent experimental measurements of 90° domain walls in PbTiOs using AFM 
by Shilo, et al. [2004]. They observed that there is a large variation in domain wall thickness with respect to 
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position (0.5 — 4.0 nm). This is not very surprising as vacancies may interact differently with different types 
of domain wahs. In our calculations, we observe that oxygen vacancies have a significant effect on the detailed 
structure of 180° domain walls. Our result is in agreement with a recent continuum study of interaction of 
oxygen vacancies with domain walls [Xiao, et al., 2005]. Xiao, et al.'s continuum model predicts that 180° and 
90° domain walls have quite different interactions with oxygen vacancies. 

There are several works in the literature on the effect of strain on ferroelectric domain walls showing that 
strain can have important effects on domain walls. For example, both experiments [Shang and Tan, 2001] and 
ab initio calculations [Shimada, et al., 2008] show that shear stress applied to 90° domain walls in PbTiOs 
develops polarization reorientation through a domain wall movement perpendicular to itself known as stress- 
induced domain switching. To our best knowledge, there are no atomistic calculations in the literature on the 
effect of strain on the structure of 180° domain walls in PbTiOa. 

This paper is structured as follows. In §2 we briefly review the previous studies of ferroelectric domain walls. 
We then present the main ideas of anharmonic lattice statics analysis of perfect and defective ferroelectric 
domain walls under strain. In §3 we report some numerical results for 180° domain walls in PbTiOs using a 
shell potential. Conclusions are given in §4. 

2 Ferroelectric Domain Walls 

Ferroelectric domain walls have been studied extensively both theoretically and experimentally. For recent ex- 
perimental investigations see [Shilo, et al., 2004; Franck, et al., 2006] and references therein. Shilo, et al. [2004] 
studied the structure of 90° domain walls in PbTiOa using atomic force microscopy (AFM). They measured 
the surface topography of a given sample and compared it with a displacement field that is obtained from the 
Devonshire-Ginzburg-Landau phenomenological model. They calculated the thickness of two domain walls at 
two different positions and observed thicknesses of 0.5 nm and 4.0 nm for the two domain walls. They then con- 
jectured that presence of point defects is responsible for this variation. Lee, et al. [2005] used a two-dimensional 
square lattice model that has a continuum order parameter interacting with a lattice of Ising spins, where the 
Ising spins model the presence or absence of point defects. They show that depending on the parameters used 
in their model, one can reproduce the main features of the experimentally observed variation in domain wall 
thickness in Shilo, et al. [2004] 's experiments. However, one should note that Shilo, et al.'s results may not 
apply to all types of domain walls as they studied only the 90° domain walls. 

Ab initio calculations of He and Vanderbilt [2001] show that oxygen vacancies have a tendency to move to do- 
main walls and pin them. They also observed that defective domain walls are Ti-centered. In their calculations, 
they had to assume a periodic array of domain walls with a high density of charge-neutral oxygen vacancies. 
There have also been other ab initio calculations of point defects in PbTiOs in the bulk [Padilla, et al., 1996; 
Park and Chadi, 1998; Cockayne and Burton, 2004]. Xiao, et al. [2005] studied the effect of oxygen vacancies 
on the structure of domain walls in tetragonal BaTiOa using a continuum theory that takes into account the 
fact that fcrroelcctrics are wide-band-gap semiconductors. In their numerical calculations they observed that 
180° and 90° domain walls behave differently in response to oxygen vacancies. In particular, they saw charge 
accumulation near 90° domain walls with a potential drop across the wall while these were absent in the case 
of 180° domain walls. 

Anharmonic Lattice Statics of Domain Walls Method of lattice statics was introduced by Matsubara 
[1952] and Kanazaki [1957] and was extensively used by Born and his co-workers [Born and Huang, 1988]. For 
more details and history see Ortiz and Phillips [1999]; Yavari, et al. [2007a,b]. In this paper, we study the 
possibility of domain wall broadening by oxygen vacancies in the case of 180° domain walls in PbTiOs. We also 
study the effect of shear and normal strains on both defect-free domain walls and domain walls with oxygen 
vacancies. 

We present a semi-analytical solution of structure of 180° domain walls using an anharmonic lattice statics 
method [Yavari, et al., 2007a,b]. From ab initio calculations of oxygen vacancies in PbTiOa [He and Vanderbilt, 
2001] and also molecular dynamics simulations of CaTiOa [Calleja, et al., 2003], we know that oxygen vacancies 
have a tendency to move to the domain walls and pin them. Therefore, we study the structure of domain walls 
with oxygen vacancies sitting on the wall. To be able to solve the discrete governing equations analytically we 
need to assume some periodicity for the collection of vacancies on the wall. We consider a 180° domain wall in an 
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infinite crystal and assume that oxygen vacancies are periodically arranged on the domain wall. Thus, the only 
restrictive assumption is the periodicity and high density of vacancies on the domain wall. In reality, oxygen 
vacancies can be distributed randomly and have a smaller density. However, the results from our calculations 
can still provide some important quantitative information on the effect of oxygen vacancies on the structure of 
180° domain walls in the nanoscale. 



Defect-Pree 180° Domain Walls Wc use a shell potential for modeling PbTiOa [Asthagiri, Wu, Choudhury, and Cohen, 
2006] and all the calculations arc performed for T = K. Denoting the collection of cores and shells by £, 
i e £ represents a core (or shell) in the collection. In a shell potential, total energy has the following form 
[Dick and Overhauser. 1964; Sepliarsky and Cohen, 2000; Sepliarsky, et al., 2004] 

^ {W}.ec) = ^^hort ({x»},,^) + fiong ({x^},e J + hen ({x^},,^) , (2.1) 

where {x'}^g^ is the current position of cores and shells. Short range energy depends explicitly on the position 
vectors of the massless shells. Long range energy is the Coulombic energy of all cores and shells, excluding 
core-shell interaction in the same atom. The core-shell energy prevents the shell collapse in each atom and is 
usually a polynomial function of the pairwise distance of core and shell in a given atom. In the equilibrium 
configuration B = {x'}jg£ C MP, energy attains a local minimum, i.e. 

8F 

-- = V i e £. (2.2) 
ax' 

However, in the case of a defective crystal, the problem is that we do not know the relaxed configuration a priori. 
Thus, we start with a reference configuration Bq = {xo}jg£ that is not necessarily force free [Yavari, et al., 
2007a]. In the case of a 180° domain wall, Bq is a nominal defect, i.e., a configuration in which cores and shells 
on the left and right sides of the wall have their positions in the bulk configurations corresponding to Ps and 
—Ps, respectively, where Ps is the spontaneous polarization (see Fig. 2.1). This reference configuration is not 
force free. Let us denote the discrete field of unbalanced forces by 



(Bo) } . (2.3) 



d: 



X' 



In anharmonic lattice statics one finds the discrete deformation mapping if : Bq ^ B that takes the chosen 
nominal defect to its relaxed configuration. Note that a different choice of reference configuration Bq leads to a 
different discrete deformation mapping (p' : Bq ^ B. If the two reference configurations are 'close', i.e., they are 
in the same energy well, we would converge to the same relaxed configuration. The discrete deformation map 
is found semi-analytically as follows. 

Taylor expanding (2.2) about the reference configuration and ignoring terms higher than quadratic in dis- 
placements, one obtains the following linearized governing equations 

or 

j 

where — 'xP — Xq. For a defect- free domain wall, atoms (cores, shells) of the same type parallel to the 
wall will have the same displacement vectors. This symmetry simplifies the linear equations (2.4) considerably 
[Kavianpour and Yavari, 2009]. Assume that the defective crystal C has a 1-D symmetry reduction, i.e., it can 
be partitioned into two-dimensional equivalence classes, i.e. 

M 



C=[}[}SIo.^^). (2.6) 



aeZ 1=1 



2 Ferroelectric Domain Walls 4 



V 



; Q o Q 


/ — domain wall 


• • • : • 1 

.0 1^, 




• • • • 


• • • t 


: ^ : 


° . ° . ° ■ ° + 
• • • 1 


• • • 
■ ■ ■ 


O O O 


T * * * 


pt ; 


i-p. 




Pt ; 


i-p. 



• Pb PbO Plane 
o 01 



■ Ti 

o 02 TiO^ Plane 
• 03 



axial strain (e ) 



A 



shear strain (v ) 




infinite crystal 



\ s ' 



2Na 



Figure 2.1: Reference configuration for a Ti-centered 180° domain wall in PbO and Ti02 planes. Note that cores and shells on 
the domain wall have no relative shifts and cores and shells on the left and right sides of the wall have opposite relative shifts, a 
and N are lattice parameter and size of the simulation box, respectively. 



where Sia{i) is the equivalence class of all the atoms of type / and index a with respect to atom i. Here we 
assume that £ is a multilattice of AI simple lattices (for PbTiOa, M = 10). Particle i is an arbitrary core (shell) 
and assuming that it is in the nth equivalence class of its type, S]a(i) is the set of those cores (shells) of type / 
that are in the equivalence class n + a. For a free surface, for example, each equivalence class is a set of cores 
(shells) lying on a plane parallel to the free surface. Using this partitioning one can write 



M 



J a = -oo 1=1 j£Sia,(i) 



(2.7) 



where the prime on the second sum means that the term a = 0, 1 
governing equations are written as 



is omitted. The linearized discrete 
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(2.8) 



(2.9) 



Unit cell displacement vectors are defined as X„ = (u^^, . . . ,u*''^)^. Now the governing equations in terms of 
unit cells displacements arc 



(2.10) 
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where A„(n) g ]g3A/x3A/^ X„,F„ G M.'^^^ and to is the range of interaction of unit cells and m = 1 would 
be accurate enough for the shell potential [Yavari, et al., 2007b]. Eq.(2.10) is a linear vector-valued ordinary 
difference equation of order 2m with variable coefficient matrices and the unit cell force vectors and the unit 
cell stiffness matrices are defined as 



( Knc 
K21C 



Ki2q 
K22Q 



K2MQ 



A/1q 



K 



M2a 



K 



n G Z. 



(2.11) 



MM 



Note that (n) explicitly depends on n and this reflects the fact that close to the domain wall force constants 
may change. 



180° Domain Walls With Oxygen Vacancies It is known for quite sometime that for a large enough 
number of oxygen vacancies in perovskites one would see a self-organized planar arrangement of oxygen va- 
cancies (see Scott and Dawber [2000]; Gopalan, et al. [2007] and references therein. See also Zhang [2004] for 
atomistic calculations and discussions on oxygen vacancies in BaTiOs and their different possible arrangements). 
Therefore, similar to the ab initio calculations of He and Vanderbilt [2001], we assume that oxygen vacancies 
interact with a domain wall and have a planar structure coinciding with the domain wall. 

Eq.(2.10) is the governing equation for a defective domain wall as long as the oxygen vacancies are arranged 
periodically on the domain wall. However, Aa{n) and F„ would change, in general. We can simplify the solution 
of the discrete boundary- value problem even further by noting that the displacements of cores and shells on the 
left side of the wall are related to those on the right side of the wall, i.e., the wall is a reflection symmetry plane. 
This reduces (2.10) to an ordinary difference equation on N2, i.e., 

A_i X„_i + Ao X„ + Al X„+i = F„ n > 2, (2.12) 

with the following boundary equations 

A_i(0) X_i+Ai(0) Xi =Fo, (2.13) 

Ao(l) Xi+Ai(l) X2 = Fi. (2.14) 

Note that Aa{n) = Aa, n > 2 and also X_„ — — X„ Vn > and hence Xq = 0. Note also that because 
of symmetry the z-componcnt of all displacements are zero and hence X„ G ]R^°, |7i| > 1. For a defective 
Ti-centered domain wall with 01 or 03 atoms removed, Xq G IR'* (see Fig. 2.2). 

The governing discrete boundary-value problem is closed by requiring boundedness of displacements at 
infinity, i.e. 

lim II X„ ||< 00. (2.15) 

One should note that the infinite crystal may undergo a rigid translation after relaxation. Condition (2.15) does 
not exclude this possibility. Note also that assuming that the domain wall is a mirror plane for displacements 
removes the translation invariancc of the governing equations, i.e., those cores and shells that lie one the domain 
wall in the reference configuration remain on the domain wall after relaxation (i.e. Xq = 0). 

When there are oxygen vacancies on the domain wall, one needs to look at boundary equations carefully. 
When modeling PbTiOa by a shell potential, an oxygen vacancy means removal of the core and shell of the 
oxygen atom and because we assume a charge neutral oxygen vacancy, there will be a charge redistribution. 
This will affect the shell charges of the neighboring shells. It is known that charge redistribution is highly 
localized. Thus, in our calculations we distribute the charge AQ = Qs + Qc, where Qs and Qc are oxygen shell 
and core charges, only to the (fourteen) first nearest neighbors of each oxygen vacancy. Most of the existing 
shell potentials have fixed charges. However, in a more realistic model shell charges should be variable to be able 
to adjust themselves to the environment, e.g., in the presence of a vacancy. There have been several efforts in 
the literature on building empirical charge- variable models for different systems based on Rappc and Goddard's 
charge equilibrium method [Rappe and Goddard, 1991] and its extensions. Here, we use a fixed-charge shell 
potential for PbTiOs and are not aware of any charge-variable shell model for this material. Therefore, we 
have to model the charge transfer due to oxygen vacancies approximately. However, this should not have severe 
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Figure 2.2: Unit cell indexing in the reference configuration of a Ti-centered 180° domain wall. Note that because of symmetry 
displacements in the z-direction are all zero. 

effects on the results as our numerical tests show. We studied the sensitivity of solutions to the exact way 
of charge redistribution; we compared two cases: in the first case we distributed the charge equally to the 
nearest neighbors and in the second case we distributed the charge to nearest neighbors depending on their 
distances from the vacancy and assuming that charge distribution is exponentially decaying. We did not see 
much difference and thus in this work we distribute the charge equally between the nearest neighbors. 

We consider three types of defective domain walls: (i) 02-defective, (ii) 01-defective, and (iii) 03-defective. 
An 02-defective domain wall is Pb-centered. See Fig. 2.1 for this notation. We see that the anharmonic lattice 
statics iterations do not converge in this case, i.e., this is not a stable configuration^. This is in agreement 
with ab initio calculations [He and Vanderbilt, 2001] that predict Ti-centered defective domain walls. 01 and 
03-defective domain walls are Ti-centered with no 01 and 03 cores and shells, respectively, on the domain wall 
in the reference configuration (see Fig. 2.2). 

The linearized governing equations can be solved exactly using the method proposed in [Yavari, et al., 
2007a]. Thus, we are able to solve the linearized governing equations exactly. Now to obtain the fully nonlinear 
solutions we use a modified Newton-Raphson iteration. Solving the linearized problem, we modify the reference 
configuration by imposing the harmonic displacements and then calculate the new unbalanced forces exactly 
using the interatomic potential [Yavari, et al., 2007a]. Continuing in this manner, if there is an equilibrium 
configuration close to the chosen reference configuration, unbalanced forces converge to zero. In the present 
work, convergence means that all forces have magnitudes less than 0.05 eV/A. Now, let us briefly explain the 
modified Newton-Raphson method, called the quasi-Newton method, that we use throughout this work. 

Newton method is based on the following quadratic approximation near the current configuration B^: 

£ (S^- + j") = £ (S''-) + V£ (B^-) • + . H (B'^) • + o (l^Y) , (2.16) 

~ k ~k 

where S = B^^^ — B^ and H is the Hessian matrix. By differentiating the above formula with respect to 5 , 

^ k 

we obtain the Newton method for determining the next configuration B^^^ = B^ + 8 : 

Ik ^ ^gfe^ . i^gk^ (2.17) 

Note that in order to converge to a local minimum, the Hessian must be positive definite. 

If the calculation of the Hessian in each iteration becomes numerically inefficient (like the present problem) , 
one can use the quasi-Newton method. The main idea behind this method is to start from a positive-definite 
approximation to the inverse Hessian and to modify this approximation in each iteration using the gradient 
vector of that step. Close to the local minimum, the approximate inverse Hessian approaches the true inverse 



^This is also the case when the defective domain wall Is under strain. 
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Hessian and we would have the quadratic convergence of the Newton method [Press, et ah, 1989]. Here we use 
the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm [Press, et al., 1989] for generating the approximate 
inverse Hessian: 

~6'' ® S'' ■ A) (g, ■ A) 



= + 7i^^ AT A + • • A) " ® (2-18) 



where C' = fH') \ A = V£'+'^ - V£\ and 



Cs y - A AT . C« . A 



(2.19) 



Calculating C*+^. one then should use C*+^ instead of H ^ to update the current configuration for the next 

- k 

configuration S'^+i = B'^ + d . If C'+^ is a poor approximation, then one may need to perform a linear search 
to refine 5*^+^ before starting the next iteration [Press, et al., 1989]. 



180° Domain Walls Under Strain Wc put both the perfect and defective domain walls under normal and 
shear strains (see Fig. 2.1). We consider both compressive and tensile strains and also shear strains both along 
and opposite to the polarization directions. We apply strains to the defective lattice by imposing displacements 
in the proper directions far enough from the domain wall (displacements are in x- and y-directions for normal and 
shear strains, respectively). Note that the defective lattice is translated rigidly on both sides of the domain wall 
outside the computational box. Applying strain to the lattice should be done gradually. For both perfect and 
defective domain walls, we start with the strain-free relaxed domain wall configuration B. Then, we apply proper 
boundary displacement to increase (or decrease) strain by the value Ae to obtain the strained configuration 
Bac- Note that Ae should be small enough such that B and Sac are close to each other. In the present work, 
Ae = 0.001 for normal strains and Ae = 0.0005 for shear strains. Next, we start with Bac and calculate B2Ae- 
Repeating this procedure, one can apply large strains to the defective lattice and obtain its relaxed configuration 
Se for a given strain e. Using the notation of the previous sections, we summarize our algorithm for applying 
strain as follows. 





Input data: B^, Ae (strain increment) 


[> 


Initialization 




> Apply B.C. to Be and obtain B^ 




> Hi = H|e=ei and set = (Ri)"! 


o 


Do until convergence is achieved 




1> Calculate forces 




t> Use quasi-Newton method to calculate C^-'+i 




> Use C'^+i to obtain B''+^ 


o 


End Do 


o 


End 



3 Numerical Results 

Perfect Domain Walls In our numerical examples we show the anharmonic displacements with respect to the 
reference configuration. The first step is to calculate unbalanced forces. Our choice of reference configurations 
makes the unbalanced forces localized in the direction perpendicular to the domain wall. However, one should 
note that a domain wall is an extended defect and unbalanced forces are not localized in the tetragonal c- 
direction. We are able to handle this nonlocality issue using the symmetry reduction idea. Since core and shell 
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Figure 3.1: a- and c-structure of a perfect Ti-centered 180° domain wall under axial strain ex- Ua and Uc are displacements along 
a-direction and the tetragonal c-direction, respectively. 



displacements are close, we only report the core displacements in this section. We assume 2N unit cells in the 
simulation box (see Fig. 2.1). Our numerical experiments show that TV = 10 would be enough for calculating 
displacements as the structure does not changed by using larger values for N . 

Fig. 3.1 shows the displacements of Pb and Ti cores in the a- and c-directions [ua and Uc, respectively) 
for a perfect Ti-centered domain wall under axial strain ex- We plot the displacements for different values of 
the axial strain between —0.03 to 0.03. As expected, Ua has a larger variation than Uc under the axial strain. 
However, note that axial strain does not have a significant effect on domain wall thickness; it is seen that all 
the distortions occur within two lattice spacings on each side of the domain wall, i,e. domain wall thickness 
is about 1.0 — 1.5 nm regardless of the value of the axial strain. Fig. 3.2 shows displacements of the same 
domain wall under different values of shear strain ^xy- The c-displacements are in the polarization direction 
for negative values of ^xy and are in the opposite direction for positive values of "ixy- Hence, we do not see 
symmetric displacements with respect to the unstrained structure. Also, we observe that unlike axial strains, 
shear strains have a considerable effect on both Ua and Uc- Similar results for a perfect Pb-centered domain 
wall are shown in Figs. 3.3 and 3.4. 



Domain Walls v^rith Oxygen Vacancies We report the displacements of a Ti-centered 180° domain wall 
with oxygen vacancies under axial and shear strains. Since displacements of 01-defective and 03-defectivc walls 
are similar, here we only report the results for 01-defective walls. Fig. 3.5 shows the displacements for an 01- 
defective domain wall under normal strain. 01 vacancies lie on the domain wall and because of symmetry have 
zero displacements, i.e. they will stay on the domain wall after deformation. For the strain-free configuration, 
all distortions parallel to the domain wall (c-displacements) occur within two lattice spacings on each side of 
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Figure 3.2: a- and c-structure of a perfect Ti-centered 180° domain wall under shear strain "fxy Ua and Uc are displacements 
along a-direction and the tetragonal c-direction, respectively. 

the wall, i.e., thickness of the domain wall in c-direction is not affected by oxygen vacancies. However, it is seen 
that structure is significantly different from that of a perfect domain wall. It is also seen that unlike perfect 
domain walls, a-displacements have the same order of magnitude as the corresponding c-displacements. We 
observe that the a-displacements are nonzero within three lattice spacings on each side of the wall. Thus, the 
thickness of an 01-defective domain wall is about 1.5 — 2.0 nm, i.e. oxygen vacancies increase the domain wall 
thickness by about fifty percent. Here similar to perfect domain walls, we see that normal strains do not have 
a significant effect on the thickness of the domain wall. Fig. 3.6 shows the displacements of an 01-defective 
domain wall under shear strain. Again, we see that shear strain has a significant effect on the displacements of 
the domain wall. Also note that domain wall thickness increases in a-direction, but this increase is less than 
the increase observed for perfect wall (see Fig. 3.2). 

The only restriction in our calculations is the high density of oxygen vacancies (similar to the existing ab 
initio calulations). The resulting stiffness coefficient matrices become highly ill-conditioned by increasing the 
period of vacancies. Therefore, we did the calculations only for the three cases where all 01, 02, or 03 oxygen 
atoms are removed from the domain wall. In all the three cases unbalanced forces in the tetragonal c-direction 
are nonzero only in two unit cells on each side of the wall. We checked this for several lower density arrangements 
of oxygen vacancies and observed the same local property for unbalanced forces. In the case of perfect domain 
walls, unbalanced forces perpendicular to the wall are very small and nonzero only in two layers on each side 
of the wall. In the case of defective domain walls, these unbalanced forces are nonzero in three unit cells on 
each side of the wall. This was also the case for several other lower density arrangements of oxygen vacancies 
on the wall. Thus, we believe that the high density of oxygen vacancies does not have a significant effect on the 
thickness of the defective domain wall, although it affects the structure. In other words, lowering the density of 




0.02 



0.01 



a 



-0.01 



-0.02 



0000000 



OOOOOOOO e ^ 



Ti 




oooooooo 

■ X W W K W W W W 



0000000 



-10 



-6 -4 -2 2 4 6 8 10 

n 



0.4 



0.2 



-0.2 



-0.4 

















Ti 







-10 -8 -6 -4 -2 2 4 6 8 10 

n 



Figure 3.3: a- and c-structure of a perfect Pb-centered 180° domain wall under axial strain e^. Ua and Uc are displacements along 
a-direction and the tetragonal c-direction, respectively. 



oxygen vacancies we expect to see changes in structure but no significant change in thickness. 



4 Concluding Remarks 

In this paper, we presented a semi-analytic study of the effect of oxygen vacancies and strain on the structure of 
180° domain walls in PbTiOs using a shell potential. We considered both Pb-centered and Ti-centered domain 
walls with oxygen vacancies sitting on them. We observed that Pb-ccntcrcd domain walls with oxygen vacancies 
are not stable (even under strain) and this is in agreement with the recent ab initio calculations that predict Ti- 
centered defective domain walls. To be able to solve for the structure of defective domain walls semi-analytically 
we have to work with a high density of oxygen vacancies on the domain wall (similar to the existing ab initio 
calculations). However, we believe that density of oxygen vacancies does not have a noticeable effect on the 
thickness of the defective domain wall. 

We observe that oxygen vacancies change the structure of a 180° domain wall significantly. One important 
effect of oxygen vacancies is that in the presence of oxygen vacancies, displacements perpendicular to the domain 
wall are of the same order of magnitude as the displacements in the tetragonal c-direction. In the anharmonic 
lattice statics iterations we observe that a-displacements have a fairly long tail, about five lattice spacings on 
each side of the wall. However, the thickness of a defective 180° domain wall is about 1.5 times that of a perfect 
domain wall. This is different from the results of a recent experimental study of 90° domain walls in PbTiOa 
using AFM and the observed large variations of domain wall thickness in the internal 0.5 — 4.0 nm [Shilo, et al., 
2004]. However, our results are in agreement with those obtained from a recent continuum model that predicts 
different behaviors of 180° and 90° in response to point defects [Xiao, et al., 2005]. We studied the effect of 
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Figure 3.4: a- and c-structure of a perfect Pb-centered 180° domain wall under shear strain "fxy Ua and Uc are displacements 
along a-direction and the tetragonal c-direction, respectively. 

strain on both perfect and defective 180° domain walls. We observed that normal strains have a greater effect on 
a-displacements but shear strains have a significant effect on both a- and c-displacements. Finally, we observe 
that the domain wall thickness does not change significantly under normal or shear strains. 
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